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^ Abstract 

m 

We define the energy- momentum tensor on lattice lor the Xcb 4 and lor the 
00 . . . ..... 

nonlinear u-model Symanzik tree-improved actions, using Ward identities or 

an explicit matching procedure. 

The resulting operators give the correct one loop scale anomaly, and in 

the case of the sigma model they can have applications in Monte Carlo simu- 

lations. 
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I. INTRODUCTION. 



The lattice regularization is the preferred framework for the non perturbative study of a 
field theory by means of numerical simulations. It gives the bonus of explicit preservation of 
gauge invariance, which is believed to be the key to understand very important phenomena, 
as for example confinement. The lattice however breaks other symmetries, as the translation 
(Poincare) invariance, and it is important to verify the correct restoration of these properties 
in the continuum limit. 

A classical symmetry is expressed at quantum level by a series of Ward identities involving 
the associated current operator. These identities can't be satisfied in presence of a non 
invariant regulator A. If the invariance must be preserved, we can expect that, after a 
suitable 0(h) operator redefinition, these relations will be broken only by "irrelevant" (that 
is to say, which vanish removing the regularization) terms. 

In general this program cannot be accomplished, as a classical invariance may not be 
compatible with the quantization. A well known example is the dilatation invariance, that 
independently of the regularization chosen cannot be restored in the renormalized theory. 
However, if an invariant regulator can be found, we must be able to restore the Ward 
identities in another scheme, whatever it would be. 

This is the case of the translation invariance, for which an invariant scheme exists (for 
example the dimensional regularization). The relative Ward identities are relations between 
insertions of the energy-momentum tensor (EMT) that can be written as 



In a concrete lattice simulation the cutoff A = 1/a is obviously finite. This gives rise to 
systematic errors in the estimation of observable quantities that can be reduced using very 
large lattices. An alternative possibility, proposed by Symanzik [8] is to redefine the lattice 
action: we can exploit the ambiguity inherent to the lattice transcription and perturbatively 
eliminate the cutoff effects to a given power of a. 

We have explicitly calculated the finite corrections necessary for the one loop definition 
of the EMT with a Symanzik tree improved action. First, as an academic exercise, we 
have analysed the simple A</> 4 model [9]. Next we have considered the nonlinear 0(N) 
bidimensional a model [10]. This is of much interest as it shares some important properties 
with lattice gauge theories [11], first of all the asymptotic freedom which is an argument to 
justify a perturbative approach. 

Since the EMT is related with the dilatation current, which as we said is broken at 
quantum level, we have also verified that the so defined operator gives the correct one loop 
anomaly. 



n 






2 



II. SYMANZIK IMPROVEMENT. 



We use the A</> 4 model to explain in detail the procedure. Starting from the classical 
continuum action S we can write a lattice regularized version of it Si att in infinite ways, 
requiring only that 

]im Si att [<f>] = S[<f>]. (2) 

a^O 

The most simple option is the substitution of the integral with a sum over lattice sites n, 
and of the derivatives with the finite differences A+</>(n) = [<f>(n + afi) — <p(n)]/a } 



^°« 0) M = « 4 E 



S[<P] + 0{a 2 ). (3) 



If we limit ourselves to a tree level (classical) improvement we can easily do better. Diffe- 
rences proportional to powers of a between the two actions arise from the derivative term, 
and can be compensated adding to the lattice lagrangian suitable irrelevant operators, or 
equivalently using an improved finite difference formally defined by 

d, = - a log (f + aA+) (4) 

with the second term expanded at the desired order. The classical action with the 0(a 2 ) 
artifacts removed is given by 

S}lS\<f>] = SlTtt] + « 4 E ^« 2 A:A^(n)A+A^(n) = S[<f>] + 0(a% (5) 

as can be easily verified expanding the finite differences. 

In the quantized theory it is necessary to parametrize the divergences of the continuum 
theory, for instance using dimensional regularization, and to define an improvement criterion. 
We require that the lattice one particle irreducible functions Y generated by S^/J coincide 
at small external momenta with the continuum ones, apart from o(Ti ) + o(a 2n ) corrections. 
Symanzik has shown that it is possible to write 

n I N k / \ 

S^M = S^y>\ + EE^E/^) [a D Y,a 2k Of k+D \n)) , (6) 

k=0 p=l j \ n J 

where the are k-dimensional lattice operators, and the coefficients f JP are calculable 
imposing the improvement criterion with a p loops calculation. 



A. Energy momentum tensor. 



From now on we work in the tree improved theory defined by 5^ t '°\ which is equivalent 
to an effective continuum (dimensionally regularized) action that can be written as 



eff 



d D x 



(7) 
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The Z constant can be determined imposing the equivalence of superficially divergent T 
functions. For the two points 1PI function we have on lattice (see app. A for the notations) 



rS t (o,o) = m 2 + | 



d A l g 



Zn m 
-4 + 



(4tt) 



log a m - Fqooo - 1 + 1e ~ T 



and in the continuum 



,(2) /n n ^ 2 , 9V 



4-D 



d u l 



(2tt 



D CO nt(l) = m 2 + 



m 



(47T) 



1 m 2 

- - 1 + IE + log 

e 47r/i z 



(9) 



Now we must subtract the pole part in (9) (we use minimal subtraction in the continuum), 
and dispose of the quadratically divergent piece in (8) by imposing the vanishing of self 
energy on lattice for p = m = with a suitable counterterm 



AS, 



(i-o) 

latt 



(10) 



Comparing the two results we obtain (we set the mass scale as // = 1/a for simplicity) 



Z,„ 



1 + 



2(4tt) 



[Iog47r-To-F ( 



0000 



The procedure for the determination of Z g is analogous. For the four points functions we 
obtain on lattice 



rj o 4 i(0, 0, 0, 0) = g - \g 2 j Aa«(0 = 9 



3 g 2 



2 (4tt) 



\~1e + -Poooo + ?o - log a 2 m 2 ] 



r m 



and in the continuum 



1^1(0,0,0,0) =9-7/ 



d u l 



(2tt 



D 2 cont (l)=9 



3 g 2 



2 (4tt) 



1 m 2 

-~1E- log 



so that 



Z n 



3 g 2 



2 (4tt) 



[lo g 47r-r -F ( 



0000 



(13) 



(14) 



It is now possible to define the EMT in a very simple way . In fact it can be shown [12,13] 
that, given the effective lagrangian L e jj defined by (7), it is sufficient to write a lattice 
operator which reduces in the leading order in a to 



T e JJ = Z^d^d^] - 6^N[L eff ] 



(15) 



where N[0] is an operator corrected accordingly to the prescription of the effective lagran- 
gian. We note that as we consider only the integrated T Mi/ , there are no problems of ambiguity 



4 



in the definition of the classical operator, so that we can neglect the so-called "improved 
Coleman term" [17,16]. 

We can also write at one loop order 



6 



Z 2 



z. 



- J2 a^a^ ^> + + y ™ V + ^W 4 



(16) 



where we have used the symmetric finite difference A^(ra) = (<f>(n + afi) — <p(n — ajl))/2a } 
and we have imposed the Ward identities to determine the Z finite corrections. We have 
explicitly checked that the two methods agree each other. 

The expression (16) reduces to the classical improved quantity if Z\ = Z 2 = Z4 = 
Z 5 = l,Zs = 0, and as we are considering only tree improvement we can neglect 0(a 2 g) 
corrections. The resulting Feynman rules give two and four lines vertices, T^ 2 ) and T^ 4 ). The 
relevant identities defined by (1) involve the insertion of these vertices in the two and four 
points irreducible functions (see figure (1)). We start from the two points identity, written 
in momentum space for the integrated T^: 



ra(o : p.- P )=(^ + ^-vW). 



2 dp v 2 dp^ 



;i7) 



We have already calculated the one loop correction to I^ 2 ) (see eq. (8)), while for the insertion 
we obtain 



lf>;p, -p) = 2T$( Pl -p)- 9 -j 7^4 Aa«(0 K + 2T(?(Z, -l)D latt (l) 



2 J (2tt / 

and substituing in the (17) we find the conditions 



'Z x - 1) = (Z 2 -l) = Z 3 = 0(g 2 

Z 4 = 1 + 9 - [ -——-D 2 l tt (l) 
4 J 2tt) 4 



I 2 a 2 I 4 



1 + 



To - T 2 + 2tt 2 Z { 



0000 



1 - l -T t 

3 



+ 



Ti T 3 - 

— + — 

4 12 



:i8) 



(19) 



(20) 



2(4tt) 2 

The determination of the Z 5 correction can be done considering the four points identity 

(o; 0,0, 0,0) = -vr (4) (o, 0,0,0). (21) 

We calculate the relevant insertion which, taking into account the relations (20), is of the 
form 



rW(0;0,0,0,0) = -V 5 Z 5 + 3 5 2 ^ 



d 4 l 
(27) 



+3/ 



d 4 l 
(27) 



4 Aa«(0 



a 2 ^ 



2LI + - Lit + ha - i 2 + -J 4 + m 



(22) 
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Using this result and the eq. (12) we obtain 

3 r d 4 l 



Z 5 = l + ~g 
= 1 + 3 



(2tt 



1-Aa«(0 P + 2m 2 + y/ 4 



9 



2(4tt) 



Tq — T2 + 27T 2 ZooOO 



1 - !ti 

3 



(23) 



B. The dilatation current anomaly. 

From the Ward identity for the Noether current connected to scale invariance 

4 D) = ^ (24) 
one can show that the scaling equation is related to the EMT trace as follows 

[ a ta - m ^- N ) r(JV)(?1 ' ' ' ' > PN) * r K ( ° ;Pl ' ' ' ' > PN) (25) 

Now, the mass derivative is equivalent to the insertion in T of a renormalized mass operator, 
and analogously the field number N can be generated inserting the renormalized operator 
</>||-. It follows that scale anomaly is equivalent to the insertion of the evanescent operator 

« 4 E(^ + ™¥ + ^). (26) 

Using the defined EMT we obtain the same results one finds with the not improved action [4]. 
This is correct, as the anomaly is proportional to the f3 functions of the theory that, as well 
known, does not depend to this order from the regularization scheme. 



III. RESULTS FOR THE IMPROVED NONLINEAR SIGMA MODEL. 



The nonlinear O(N) sigma model is defined by the path integral 

Z = J V$(x)6(P(x) - 1) exp{-l J d 2 ydj(y) ■ dj(y) 



(27) 



where <f> = (</>i, • • • , cf>^) are N scalar fields. In the continuum we use as usual the dimensional 
regularization scheme, with minimal subtraction at the scale // = 1/a, and in order to 
eliminate infrared divergences, which are an artefact of perturbative approximation, we add 
to the action a mass term which breaks the 0(N) symmetry. With the chosen regularization 
there are no contributions from the path integral measure, and we can write the action as 



s 



1 



1 



—dfi ■ dfi + ^;<V<V - —a 



m 
t 



(28) 



with the parametrization 
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7T = ^i, • • • ,<j>N-l, 

a = cf> N = v / i ~ " 



7T • 7T. 



(29) 



This lagrangian generates an infinite number of vertices, but at a given perturbative order 
only a restricted set of them gives a contribution. At one loop we need the propagator of 7r 8 - 
fields Dij(p) and the four points vertex Kjfc/(pi,P2,P3,P4)- 

The 0(N) invariance gives strong constraints to the renormalized action. Solving expli- 
citly the relative Ward identities [f 4] one can prove that this can be written in terms of two 
renormalization constant Z } Z t} 



SW = j d D xp D ~ 2 



7T Z 



On lattice the tree improved action can be written as 

^latt ^ = ^0 ~\~ S\ ~\- S meas 



(30) 



(31) 



wi 



ith 



S = a 2 Y^ 

n 



^A+vr(n) • A+7?(n) + ^A+a(n)A+ ( r(n) 



m 



t 



-a(n) 



2 2 

— A+A+7?(n) • A+A+vr(n) + — A+A>(n)A+A>(nl 

2^_f /i/iV/ /i/iV/ 1 a* /i v / /i /i v > 



1 

u vneas 



^ — log(j(n) 



(32) 

(33) 
(34) 



the last term coming from the exponentiated path integral measure. This action is equivalent 
to an effective continuum one, that must be of the form (30). To determine the finite 
renormalizations Z, Z t it is sufficient to compare the two point irreducible function on the 
continuum 



r l 2 ,LtCp, -p) = 8 u 



p 2 + m 2 



1 ( 2 N - 1 2 

— V H rn 

4tt V 2 



log 



am 
Air 



+ IE 



and on the lattice 

^\j]latt{Pi ~P) = ~Y 



P 



t 

+ 2 



d 2 l 
(2vr) 2 

,2 



Aa«(0 + 



D latt (i) £ ((p+if p + (p^i) p + ^(p+i) 4 P + 



(35) 



(36) 



Using the lattice integrals properties and suppressing irrelevant terms we can write 



rgL(p,-p)^[p 2 + - 2 + 



+ 



d 2 l 
f2vrV 



Dl att (l) 



' N ~ 1 2 / 7 2 7 sP* 

— - — m + (f + 4 cos a/i — 2 cos al\) — 



(37) 
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where the quadratic divergence in the integral has been canceled by the measure term. The 
final result is 

Z, = 1 + ±(N -2) (; E + E + log |) + ± (| - i E 5 + |£i + |£i) (38) 

Z = 1 + — (TV - 1) [ 1E + E + log — . (39) 

47T \ 47T/ 



A. Energy momentum tensor. 



At first we note that we expect to obtain a 0(N) invariant result, except for contributions 
proportional to the symmetry breaking mass or to the motion equation 



)7T 



A+A-7? 

/i /i 



7T 



(7 



>7T 



m 



a 



So at the order we are working the operator can be written as 



rp _ T-i(O) , T-i(l) , rp( m .e.) , rp(l) 
- 1 - W> - 1 - iu/ ' - 1 - in/ ' - 1 - in/ ' - 1 - in/ 



Wl 



ith 



t(o) 

/J/V 



rn[ra.e.) 

/J/V 



% _ 



-E4 



7 

• m a 



CL 

^ SSl a tt 
^47T • ^ , 



+ <V 



(40) 
(41) 

(42) 
(43) 

(44) 
(45) 



When Z\ = Z 2 = 1 the term Tffi reduces to the naive lattice transcription of the classical 
operator. T^) cancels the 0(a 2 ) lattice artifacts at tree level, and we have not considered 
0(t) corrections for it as they would be irrelevant to the order we are working. T^' 6 ^ is a 
possible contribution proportional to the motion equation, and as we said it doesn't need to 
be O(N) invariant. It is not present at tree level, so Z 4 = 0(t). T^J is an operator with only 
the hypercubic lattice symmetry that can be required to ripristinate the 0(a 2 ) rotational 
invariance. 

As in the A$ 4 model it is possible to calculate the Z corrections by using two equivalent 
methods. We start evaluating the insertion with external momenta p = k and q = (see 
the graphs 10,11,12 in fig. (1)). In this regime the Ward identity reduces to 



u1^(M,o) = -MlM°)- 



The contribution of the various insertions can be written as 



(46) 



8 



I o (k,0) = %Vt \z 4 (k 2 + 2m 



Z, n 



m 



h(k,0) 



6 tl f d 2 l 



(2*) 



2D latt (i) + (N-1) 



Dlatt(l)Dlatt(l + *Q 
D la „(k) 



D latt {l) 



A^{k + l,l)~ V [\ A Pp( k + h k) 



x [A^(k + /, /) + 



, ^ + 1 2 
H m 



A^{p, q) = v^ v + y [p^q v + p^q v ) + [n^v\ 

Summing the three contributions and using (46) we obtain 



Z< = -(*-!) 
Z„, = l-^-(A ; -l). 

07T 



(47) 

(48) 
(49) 
(50) 

(51) 
(52) 



The other constants can be calculated from the insertions of the integrated operator. With 
the same notation we have 



hip, ~P) = Y 



2Z iP^Pv + J (p»pI + pIpJ) ~ V ( Z 2p* + + Z 3 pl + Z m m 2 + 



+ 2Z 4 (p 2 + m 2 ))] 

TTl 2 f (Pi 

h(p, ~p) = f>ijti»v-r- 



h(p, ~P) = 



2 J (2tt) 2 
d 2 l 



(2*) 



:D latt (l) 



Dl tt (l)\2D latt (p + l) + (N-l)m 2 

6 



2A, 1J {p + l, p + l)--2L [A pp (l + p,l + p ) + (N + IK 



(53) 
(54) 

(55) 



Taking all together, and using the Ward identity (17) we obtain 



Z 1 = l + —(4^5 + E - 3Ex - 2E 3 - 3E 2 - 20) 

1Z7T 

Z 2 = l + t^-(4E 5 - 2E - 2E 3 + 8E 4 - 3(iV - 1) - 5) 

1Z7T 

Z 3 = -^(3£o - 3E X - 3£ 2 - 8E 4 - 12). 

D7T 



(56) 
(57) 
(58) 



This result can be checked by an explicit determination of a lattice transcription of the 
effective EMT 



T e JJ = ~N[dJ-d v f\ - S^N[L eff ], 



(59) 



where L e jj is defined by (30). It is sufficient to match lattice and continuum insertion of 
operators in the two points function. 

For reference we list the results for the matching coefficients 
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h = l + 



b 2 
b 3 



c 2 = c 3 



defined by 



47T 



32 

-4 - E 2 - E x - log h 7-B 

47T 



8tt 
t 



(jV-2)(-log£- 7B + £ 6 + l)-|[ 



10 [12-3^0 + 3^ + 3^ + 8^4 

127T 

t r 32 

6 4 = 6 5 = -l(iV - 1) ' " ■ ' 



8tt 



1 - E 6 + log h 7-B 

47T 

(TV - 2)E 6 - (N - 1) (log ^ + 7s) " " E + 



47T 



(iV-1) 



32 

log V^e - E 6 

47T 



(60) 

(61) 

(62) 

(63) 

(64) 

(65) 
(66) 



N[d^ ■ d v <f>] = WAJ ■ A v <f> - — [A J • A 3 J + A^ • A v <f>) + b 2 S^A p <f> ■ A p <f> + 

+ b^S^A^cj) ■ A^ + btS^rmfa + hS^n^^l, (67) 

a 2 8 S 

iV[«9 p ^ • d^] = Cl A p <^ • A p <^ - - A J ■ A 3 p (f + c 2 mV + csv?^. (68) 



B. The dilatation current anomaly. 



In analogy with the A$ 4 model, we find that the scale anomaly is proportional to the 
irrelevant operator 



n ^ 

whose insertions are equivalent to those of 



-m 2 a 



« 2 E W 



at 2 (77T am 2 / 



Independently of the action (improved or standard) we find 
(iV-2) 



© = « 2 E 



2tt 



-^A^- A^ + m 2 cr 



, (iV-3) 2 (AT - 1) ^ <^ a „ 
+ - — : -mV - t- — : -ir ■ 



47T 



47T 



Stt 



and we obtain the correct one loop (scheme independent) results: 

XN-2) 



m = -t 
at) = t 

7m (t) = ~t 



2tt 

(N-l) 
2vr ' 

Air 



(69) 

(70) 

(71) 

(72) 
(73) 
(74) 
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C. Results for the standard theory. 



We report the results obtained lor the no improved theory. With the same notations 
used previously we have lor the renormalized lagrangian 



Z t = l+t 
Z = 1 + t 



(N-2) 


log 


32 


47T 


47 


(iV-1) 


log 


32 


47T 


47 



+ 



The relevant insertion ol composite operators appearing in gives 



4tt 



32 

4 + log — + IE 

47T 



t 

8tt 



b 2 
b 3 
h 

Cl = 1 + 
c 2 = c 3 = t 



3 32 
2 _ vr + (iV - 2) ( - - log — - 1E 



t 



- 1 [1 - 

7T 

(iV-1) 



8tt 

47T 

(N-l) 
Air 



3 n 32 
h log h 7b 

2 6 4tt ' 



1 / 32 

_ (iV _ 2 ) - (TV - 1) (log — + 1E ) + 3tt - 6 

\ 32 1 
log h 7b 

. 4tt ' 2, 



The renormalization constants lor the energy-momentum tensor are 

Z 1 = l--L[4 + 7r], 

47T 

Z 2 = 1 + [3tt - 8 - (TV - 2)] , 

47T 

Z 3 = ~t [1 - 7T] , 
7T 

^4 = ^ [iV - 1] , 

Z m = 1 - [N - 1] . 

07T 



IV. CONCLUSIONS. 



We have defined on the lattice, to one loop, the energy-momentum tensor lor two Syman- 
zik tree improved actions, the A$ 4 scalar theory and the bidimensional nonlinear a model. 
In the a model case, due to the asymptotic freedom ol the theory, the calculation is reliable 
in the asymptotic scaling regime. 

We have found the correct scale anomaly, which as expected is independent ol the regu- 
larization scheme. 
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For the a model the corrected EMT we have found can be a starting point for a successive 
non-perturbative determination with numerical method [18]. As the T 00 component of EMT 
is the energy density, it can be used on lattice for the determination of the mass spectrum. 
We think another interesting possibility is the test of EMT in variational computations, in 
order to clarify what are in this framework the influences of lattice artefacts. 



APPENDIX A: NOTATIONS AND INTEGRALS. 

The lattice propagator is defined as 

Diaud) = j^— 2 , with / 2 = E 4 si* 2 (f • (Al) 

We use also 

^ = ~ sina/^. (A2) 

The lattice integral needed for calculations with the standard (not improved) action for A$ 4 
can be evaluated in terms of integrals of modified Bessel functions [15]. All the relevant 
cases can be expressed in term of two parameters 

h = I ,f* 1 2 =^"7^(1 + Foooo -IE- log aW) + 0(a 2 m 4 ) (A3) 
J (Z7r) 4 q z + m z a z (4^)^ 

h = I i^yW^y = (^ (Fooo ° " 1E ~ log a2m2) + 0{a2m2) (A4) 

f dH 1 1 1_ 2 

3 " J (2tt) 4 (<p + m 2 ) 3 " (4vr) 2 2m 2 + ( ° ' [ ' 

d 4 l q^q v _ V 



1 \ Oj^ Tfi^ 

,* 2 , 2 x 2 2 y ^oooo - - + 7— — t(27T 2 Z 000 + Je ~ 1 ~ ^0000+ (A6) 
[q l + m 1 ) 1 la 1 \_\ 8/ (47rJ 2 

+ log a 2 m 2 )] + 0(a 2 m 4 ) (A7) 
h = j ^yj^ys = i^(^ooo - 2.^0000 - 7* " l°g« 2 - 2 ) + O(aW) (A8) 



where 



i^oooo = / ze- 8z I*(2z)dz + / ze" 8 * 
Jo Ji 



8z 



(Attz) 



4.369 (A9) 



9 / >7r o? 4 o 1 

^oooo = a 2 / 77rh^ - 0-155 A10 

For the standard sigma model case all integrals are analitically calculable in term of first 
and second type elliptic functions. 
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St 

S3 

S4 



d 2 l 


1 1 




am 

y 


(2tt) 2 


p _|_ m 2 ~ 4tt 


D 32 


d 2 l 


cos all 1 


loj 


? ? 
am 

y 


(2vr) 2 


p _|_ m 2 4tj- 


D 32 


d 2 l 


cos a/i cos a/ 2 




1 , 

— log 


(2tt) 2 


I 2 -\- m 2 




d 2 l 


cos 2 a/i 1 


loj 


am 

y 


(2tt) 2 


p _|_ m 2 4t,- 


3 32 



+ 0(a 2 m 2 ' ] 



-- + 0(a 

4 v 

a 2 m 2 1 



2 2\ 

m 



32 7T 



+ f + 0(a 2 m 2 ; 

2_2\ 



1 1 

- H V U[a-m 

2 7T 



(All) 
(A12) 
(A13) 
(A14) 



The integrals which appear in calculations with the improved actions are not analitically 
calculable. For the A$ 4 model we need to extract the asymptotic a — > values of 



Mi 

M 2 



d 4 l n /A r m 2 Z — Zoooo „ , 2 4 , 

l^iatty-) = h — , A xo -^o + " ™ 



(2vr) 4 ^ (4tt) 2 



+ 0(a 2 m 4 



2 2\ 

m 



^3 = / (^J^fa«(0 = '3 + 0(a 2 , 

r dH Ti 2m 2 

M, = J j—y 4 D 2 att (l)q 2 = -1 + — (27r 2 Z 00 oo + 7£ " 1 " T 2 - F 0000 + 



2 2\ 

- m 



log a 2 m 2 ) + 0(a' 



m 



(A15) 
(A16) 
(A17) 

(A18) 



d 4 l 

W) 



To. 



1 



+ 7^7(^0000 - 27r 2 Z 00 oo -IE- log a 2 m 2 ) + 0{a 2 m 2 ) (A19) 



(4vr) 2 (4tt) 2 
T3 2m 2 



2_4^ 



(A20) 
(A21) 



M r = «7 (^-(O? 4 = ^ 4 + 0(a 2 m 2 ). 

The general strategy is to add and subtract integral resulting from the standard lattice 
formulation from the improved one. In this way we obtain a divergent piece that can 
be calculated explicitly (as is a standard integral), and a finite expression that can be 
evaluated numerically. We give the values of the relevant constants in Tab. I. For the sigma 
model we need 

d 2 l 



Si 

£2 
§3 
S4 

s 5 



(2tt) 2 

d 2 l 

(2t) 

d 2 l 



D,att(l 



(2tt) 2 

d 2 l 

(2*) 

d 2 l 



(2*) 



Dlatt{l 
^Dl att {l 



am 



cos a/i = — I — it — log 

4tt \ 6 32 

cos a/i cos a/o = — | -Ei + 4 — log 

47T \ 

cos 2 a/i = — ( E 3 + 4 - 2tt - log 

47T \ 

sin 4 all = — E4 + 0(a 2 ) 

47T 





(A22) 


) +0(a 2 ) 


(A23) 




(A24) 


^) + o,.>, 


(A25) 




(A26) 
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f d 2 l 1 
^6= -——zDi att (l) sin 4 ak cos ak cos al 2 = —E 2 + 0(a 2 ) (A27) 
J (Ztt) z 47T 

d 2 l „„ _ 1 



SV = i(2^ B - ( " = W +0( » 2) (A28) 

which have been evaluated with the same technique, in terms ol the constants reported in 
Tab. II 



14 



TABLES 



Z 0.129 j^T -0.008 

T t 0.045 j£pT 2 -0.005 

T 3 0.029 0.009 

TABLE I. The constants Zq and 



E -0.5928 0.0245 

0.0443 E 3 -0.3127 

0.1056 E 5 -0.0577 

£ 6 0.8706 £ 7 7.041 

TABLE II. The constants E { 
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